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We study a carbon nanotube quantum dot embedded into a SQUID loop in order to investigate the 
competition of strong electron correlations with proximity effect. Depending whether local pairing 
or local magnetism prevails, a superconducting quantum dot will respectively exhibit positive or 
negative supercurrent, referred to as a or tt Josephson junction. In the regime of strong Coulomb 
blockade, the to tt transition is typically controlled by a change in the discrete charge state of 
the dot, from even to odd. In contrast, at larger tunneling amplitude the Kondo effect develops 
for an odd charge (magnetic) dot in the normal state, and quenches magnetism. In this situation, 
we find that a first order to tt quantum phase transition can be triggered at fixed valence when 
superconductivity is brought in, due to the competition of the superconducting gap and the Kondo 
temperature. The SQUID geometry together with the tunability of our device allows the exploration 
of the associated phase diagram predicted by recent theories. We also report on the observation of 
anharmonic behavior of the current-phase relation in the transition regime, that we associate with 
the two different accessible superconducting states. Our results ultimately reveal the spin singlet 
nature of the Kondo ground state, which is the key process in allowing the stability of the 0-phase 
far from the mixed valence regime. 

PACS numbers: 72.15.Qm,73.21.-b,73.63.Fg,74.50+r 



I. INTRODUCTION 

Realizing a Josephson junction with a carbon nanotube 
as a weak link opened up the way to a new class of na- 
noelectronic devices combining both quantum confine- 
ment at the nanoscale and the Josephson effect In 
these junctions the critical current could be first thought 
to be maximized when a discrete electronic level of the 
quantum dot comes into resonance with the Cooper pair 
condensate of the electrodes, thus allowing an electro- 
static tuning of the magnitude of the critical current 
The Josephson effect in quantum dots is however more 
complex, because it is governed by the interplay of elec- 
tronic pairing and strong Coulomb interaction on the 
dot^^"^^. When superconductivity dominates, the super- 
conductor wave function spreads over the dot, inducing 
a BCS-singlet ground state, i.e. a standard Josephson 
junction (dubbed the 0-state in what follows) In the 
other extreme regime of large electron-electron interac- 
tions, the quantum dot enters the Coulomb blockade do- 
main, and its charge is locked to integer values, alter- 
ing the superconducting state. For an odd occupancy, 
quantum dots behave like a spin S = 1/2 magnetic im- 
purity that competes with Cooper pair formation, and 
the ground state can become a magnetic doublet. In this 
situation, dissipationless current mainly transits through 
a fourth order tunneling process reordering the spins of 
Cooper pairs, thus leading to a negative sign of the su- 
percurrent, which is referred to as the 7r-type Joseph- 
son junction^'^'^^. These two antagonist superconducting 
states, associated with a sharp sign reversal of the dissi- 
pationless current at zero temperature, can hence allow a 
first order quantum phase transition by tuning the micro- 



scopic parameters in the quantum dot. In the case of very 
strong Coulomb blockade, the O-tt transition is achieved 
by modifying the parity of the electronic charge on the 
dot (valence is easily changed using electrostatic gates), 
so that the supercurrent sign reversal occurs at the edges 
of the Coulomb diamonds. A more intriguing regime oc- 
curs for intermediate Coulomb repulsion (associated to 
moderately small values of the tunneling amplitude com- 
pared to the charging energy), in which Kondo correla- 
tions take place: in the normal state, the magnetic impu- 
rity of the odd charge state is screened through spin-flip 
cotunneling processes^^, providing a non-zero density of 
states at Fermi energy. This so-called Kondo resonance 
allows the Cooper pairs to flow normally in the super- 
conducting state, and a 0— type Josephson junction is 
therefore recovered^^'^^'^^'^^'^^. Here we explore in detail 
how superconducting transport is affected by the pres- 
ence of Kondo behavior and we finely tune the O-tt quan- 
tum phase transition in this more complex regime by 
controlling the quantum dot microscopic parameters. 



II. CHARACTERIZATION OF THE 
NANOSQUID 

A. Sample fabrication 

Here we investigate supercurrent reversal in a carbon 
nanotube Josephson junction using the nano-SQUID ge- 
ometry, which implements two Josephson junctions in 
parallel built with a unique carbon nanotube^. The 
single-wall carbon nanotubes were obtained using laser 
ablation and then dispersed in a pure dichloroethane 
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solution using low power ultrasounds. A degenerately 
doped silicon wafer with a 450nm layer of Si02 on top was 
used as backgate. A first optical lithography step pro- 
vided alignment marks in order to locate the nanotubes 
by scanning electron microscopy. The superconducting 
loops and the sidegates were fabricated using aligned 
e-beam lithography, followed by e-beam evaporation of 
Pd/Al bilayer (with respective thickness 4nm/50nm). All 
measurements were performed in a dilution refrigerator 
with a base temperature of T = 35mK, and the filtering 
stages were similar to the ones performed in Ref. 4. Sam- 
ples were current-biased, either for DC or lock-in mea- 
surements (with an AC amplitude of lOpA), in order to 
measure directly the switching current or the differen- 
tial resistance of the device. The nano-SQUID switching 
currents /gw were detected via a digital filter which mon- 
itors the estimated variance of the average DC voltage, 
see Appendix C and Ref. 31. Fig. la shows a SEM micro- 



phase shifts (0 or tt) of the two Josephson junctions, and 
(^ci,^c2) their respective critical currents. The critical 
current modulation is thus shifted by (/)o/2 between the 
— and the tt — SQUID configuration. 
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FIG. 1: Nano-SQUID characteristics. a) SEM mi- 
crograph of the measured nano-SQUID. The two sidegates 
(SG1,SG2) are colored in green, the nanotube defining two 
quantum dots is shown in red, and the superconducting leads 
are in blue, b) Map of the normal state zero bias conduc- 
tance dl/dVsd vs the two sidegate voltages at magnetic field 
B — 75mT, temperature T — 35mK and backgate votage 
Ubg = OV. Black triangles indicate the odd occupancy re- 
gions of the first quantum dot (QDl). A line cut of the linear 
conductance at fixed Vsg2 = —5.25V is also shown. 

graph of the measured nanotube SQUID with two 350nm 
long nanotube Josephson junctions (JJl and JJ2). Using 
the second quantum dot as a (possibly tunable) refer- 
ence junction, a precise control over the energy eo and 
linewidth V of the first quantum dot is accessible by tun- 
ing a pair of local sidegates and a backgate (Usgi, Usg2, 
V^Q respectively), see discussion below. Such a geome- 
try allows us to directly measure the Josephson current 
of a single junction via the magnetic field modulation of 
the SQUID switching current /gw, see Ref. 4. Indeed, the 
critical current of an asymmetric SQUID with sinusoidal 
current-phase relation (taken here for simplicity) can be 
written as: 



(j) (^1 + (^2 



'(/cl-/c2)'+4/el/c2 

V ^0 ^ / 

(1) 

where (j) is the flux modulation of the SQUID, 0o = h/2e 
is the magnetic flux quantum, ((^1,(^2) are the intrinsic 



B. Normal state transport properties 



Fig. lb presents the nano-SQUID stability diagram 
dl/dVsd at zero bias vs Vsgi and Vsg2 in the normal 
state (a perpendicular magnetic field of B = 75mT is 
applied to suppress superconductivity), at a given back- 
gate Vbg = OV. This diagram resembles a weakly tilted 
checkerboard pattern, which is typical for two parallel un- 
coupled quantum dots in the Coulomb blockade regime, 
with a weak crosstalk of about 4%. The line-cut at 
fixed VsG2 = —5.25V emphasizes the regions of high and 
low differential conductance associated with the Kondo 
ridges and Coulomb blockaded valleys respectively. One 
can indeed distinguish easily between even and odd oc- 
cupancies in each dot from the sequence of conducting 
and blocked regions: dark blue pockets denote regimes 
where both dots are blocked (in an even-even configura- 
tion of the double dot setup), green lines correspond to 
the situation of a single dot in the Kondo regime (see 
arrows) while the other remains blocked (in an even-odd 
configuration), and red spots show the case where both 
dots undergo the Kondo effect (in an odd-odd configura- 
tion)4'30. 

An operating region at a different backgate voltage 
Vbg = —0.3V is shown with greater detail on Fig. 2a. 
For VsGi between 1.70V and 1.95V, JJl has an odd oc- 
cupancy associated with a differential conductance close 
to 2e^ /h due to a well-developed Kondo effect. Further- 
more, JJ2 clearly has an even number of electrons for 
^SG2 between —4.85V and —5.15V, because of its small 
contribution to transport in this range. In order to show 
the influence of the backgate voltage VbGi we have plot- 
ted on Fig. 2b the differential conductance vs Fsgi for 
the odd occupancy region of JJl corresponding to the 
white cut on Fig. 2a, taking five different values of Vbg 
from —0.3V to —0.7V. By applying VbGi the sidegates 
experience a capacitive crosstalk of —21.5% and —17.4% 
for sidegate 1 and sidegate 2 respectively, as seen by 
the global shifts of the conductance traces. The applica- 
tion of a backgate voltage thus modifies the occupation 
number on the dot, but also the tunnel linewidth F^^'^^. 
Indeed, by varying the backgate voltage and correcting 
the sidegates voltages for crosstalk, it is possible to keep 
the local Coulomb repulsion U and the level position eo 
on the quantum dots relatively constant^^, while the hy- 
bridization F of the first quantum dot (QDl) experiences 
sizeable variations up to about 20%, as we discuss now. 
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FIG. 2: Kondo correlations and voltage/current char- 
acteristics, a) Zero bias conductance dl/dVsd vs the side- 
gate voltages VsGi and Vsg2 in the normal state (under an 
applied magnetic field B = 75mT) at Vbg = —0.3V, where 
JJl presents a well established Kondo ridge for 1.7V< Vsgi < 
1.95V and JJ2 has an even occupancy for — 5.15V< Vsg2 < 
—4.85V. b) Zero bias dl/dVsd conductance vs sidegate volt- 
age VsGi in the normal state along the white line in panel a) 
for five different backgate voltages Vbg between —0.3V and 
—0.7V. VsG2 was corrected for crosstalk in order to follow 
the white line, but Vsgi was shown as measured, c) Finite 
bias differential conductance dl/dVsd vs source-drain voltage 
V^d in the normal state, taken in the middle of the Kondo 
ridge of JJl (red star in panel a), corresponding to a level 
position eo = 0) for seven different backgate voltages Vbg be- 
tween —0.7V and —0.3V. The width of the Kondo resonance is 
modified by the backgate voltage while eo and U are kept con- 
stant, implying a variation of F. d) Typical superconducting 
voltage/current characteristics of the nano-SQUID for three 
arbitrary values of the gate voltages (Vsgi,Vsg2)- The data 
are analyzed throughout the paper by recording the switch- 
ing currents /sw obtained from such voltage/current plot, see 
Appendix C on the experimental technique that was used. 



C. Tuning the hybridization with the gates 

On Fig. 2c, the Kondo resonances taken in the middle 
of the odd occupancy region of QDl (see corresponding 
red star in Fig. 2a) are superimposed for different values 
of Vbg- The hybridization F can be extracted for dif- 
ferent values of the backgate voltage V^g from the half- 
width at half-maximum Vk of the Kondo resonance in 
the finite bias conductance. In order to extract system- 
atically Vk^ we used a Lorentzian lineshape with fixed 
background corresponding to the QD2 contribution to 
transport and to a small elastic cotunneling component 
for QDl^^. Qualitatively, we note the clear increase of Vk 
that is achieved by shifting the backgate voltage to more 
negative values, which is related to the gate-induced en- 
hancement of the hybridization F reported above. More 
precisely, in the scaling limit of the Kondo problem^^, a 
universal behavior of all physical observables is obtained 
as a function of the Kondo scale, here expressed as a 
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with U the Coulomb repulsion on the dot, F its to- 
tal hybridization to the leads, and eo its energy shift 
(taken by convention to zero in the middle of the dia- 
mond). This expression applies in the limit U ^ T (for 
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FIG. 3: Hybridization T vs backgate voltage Vbg- This 
determination is performed in the middle of the Coulomb di- 
amond of JJl, see text for details. A control with 20% ampli- 
tude variations of F is thus achieved by tuning the backgate. 
For visibilty error bars are not indicated here, but are plotted 
later on Fig. 6. 

— U/2 < eo < and contains a yet undetermined 

prefactor a, which depends on the physical quantity un- 
der consideration. To obtain the value of the charging 
energy /7, we have considered the Coulomb stability dia- 
gram of JJl, see Appendix B 1, and extrapolated the dia- 
mond edges to large bias. Because of the large linewidth 
of our strongly coupled nanostructure, the determina- 
tion of U leads to moderate error bars, and we estimate 
U = 0.80 ± 0.05 meV. Now focusing on the differential 
conductance dl/dVsd from now on, a more precise defi- 
nition of a is set by our choice of Vk as the half width 
at half maximum of the finite bias Kondo peak. In a 
near equilibrium situation (corresponding to very asym- 
metric barrier to the left and right leads) and in the 
regime U ^ 6F, we find that the unknown parameter 
is given by a^^ 2.8 from Numerical Renormalization 
Group (NRG) calculations^^. However, in our experi- 
ment the conductance is tuned to its maximum value 
Ae^ /h {i.e. 2e^/h per dot), corresponding to equally 
balanced tunneling amplitudes from each leads (we note 
that values slightly above Ae^ /h can be achieved in our 
double-dot device, which we attribute to small extra elas- 
tic contributions from symmetry-broken orbital states of 
the carbon nanotube). In that case, decoherence of the 
Kondo anomaly is induced by the antagonist pinning of 
the Kondo resonance to the split Fermi levels of each 
lead, reducing the half-width Vk of the finite-bias con- 
ductance peak compared to the equilibrium situation. 
For the relevant regime U ^ 6F, one can estimate^^ a 
reduction by 50% of the linewidth, so that we finally fix 
a = 1.4. Because there is to date no fully controlled the- 
ory of the finite bias Kondo resonance, we believe that 
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the unprecise choice of a wih introduce the largest error 
in our determination of F, and hence of the phase bound- 
ary analyzed in Sec. Ill B. The final backgate dependence 
of the hybridization T is shown in Fig. 3 for the middle- 
point (particle-hole symmetric) of the Coulomb diamond 
of JJl. The variations of F with the backgate Vbg are 
quite sizeable (up to 20%), according to the exponential 
dependence of the Kondo scale (2) and constitute a cen- 
tral piece of the analysis in the superconducting state, 
allowing to span a large part of the phase diagram of 
the O-TT transition. We also stress that changing the lo- 
cal sidegates non-only allows to tune the energy levels 
in the dots, as is clear in Fig. lb), but also modifies the 
hybridization F. The complete evolution of F with back- 
gate and sidegate voltages can be tracked by the analysis 
of the Kondo anomalies using Eq. (2), and leads to a 
greater range of variations (up to 50%). 



III. EXPERIMENTAL STUDY OF THE FIRST 
ORDER O-TT TRANSITION 

Having characterized the normal state properties of 
our device, we now focus on the superconducting behav- 
ior of the nano-SQUID. Fig. 2d shows typical voltage- 
current characteristics obtained at three arbitrary gate 
voltages in the superconducting state. For all setpoints 
that we measured, the nano-SQUID shows an abrupt 
transition to the finite voltage branch indicating an un- 
derdamped device, with an hysteretic voltage-current 
characteristics^'^. The current at which this sharp jump 
occurs defines the switching current /gw, which can be 
precisely determined via a digital filter^^ calculating the 
maximum variance of the measured dc voltage, see Ap- 
pendix C. Switching currents of approximately 3pA up 
to a few nA can thus be detected in a fully automated 
fashion. 



Comparison between valence induced and 
Kondo induced — tt transition 



Here we focus on a comparison of the O-tt transition 
behavior in two different correlation regimes, achieved 
in two distinct regions of the sidegate checkerboard dia- 
gram. Fig. 4a and Fig. 4d show both operating regions: 
a) corresponds to the situation already studied in the 
normal state, where fully developed Kondo correlations 
take place for the odd region of JJl, while d) reveals an 
odd charge state of JJl where Kondo correlations do not 
arise (Kondo temperature smaller than the base temper- 
ature of the cryostat). These distinct physical regime 
are similarly witnessed on panels b) and e), which show 
the normal state conductance trace along the white line 
in the conductance maps. In c). Coulomb blockade is 
fully overcome by the Kondo effect in the odd region of 
JJl, while in e) Coulomb blockade is robust throughout 
the entire gate range. The most interesting comparison 
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FIG. 4: — TT transition in Kondo and Coulomb block- 
aded odd-charge states. Panels a)-b)-c) respectively show 
the stability diagram, the normal state conductance along the 
white line in a), and the supercurrent, in the case of a well- 
developped Kondo effect. Panels d)-e)-f) respectively show 
the stability diagram, the normal state conductance along 
the white line in d), and the supercurrent, in the case when 
Coulomb blockade is not overcome by the Kondo effect. Su- 
percurrent sign reversal, observed by the crossing of the two 
curves at = and = 00 /2 respectively, penetrates in 
c) deep within the odd charge Coulomb diamond thanks to 
Kondo screening, in contrast to f), where the O-tt crossing 
occurs precisely when valence changes on the dot. 



between the two regimes occurs in the superconducting 
state, see panels c) and f). In panel c), a supercurrent 
reversal (indicated by the crossing of the two curves as- 
sociated to two different magnetic ffux, see Eq. 1), oc- 
curs within the odd charge state of JJl, showing that the 
Kondo effect plays a crucial role in triggering the — tt 
transition. In contrast, panel f) shows that the super- 
current changes sign concomitantly with the increase of 
valence of the dot, in agreement with expectations in the 
strong Coulomb blockade regime^^. 



B. Tuning the O-tt transition with controlled 
changes in the Kondo temperature 

As shown in Sec. II C, it is possible to tune the hy- 
bridization F with the gates, which we exploit to charac- 
terize more globally the O-tt transition phase boundary. 
Indeed, we saw previously that Kondo correlations in JJl 
are strengthened when Vbg goes from —0.3V to —0.7V 
in the operating region of Fig. 2a. In order to explore 
precisely the influence of Kondo correlations on the O-tt 
transition. Fig. 5 presents six different plots of /sw ^6^- 
sus VsGi (along the white line on Fig. 2a) at different 
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backgate voltages and magnetic fields. 
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FIG. 5: Switching current behavior of the nano- 
SQUID in the Kondo regime. The local gate voltage 
VsGi dependence of the switching current /sw is recorded for 
five different magnetic fields strengths (plotted as different 
colored lines). The various panels are displayed by following 
the decrease of the backgate voltage Vbg from -0.3V to -0.7V 
(the case Vbg = —0.5V was previously shown in Fig. 4a,b,c). 
This decrease of the backgate voltage allows to strengthen 
the Kondo effect, progressively shrinking the domain associ- 
ated to TT-junction behavior. Note that Vsg2 was corrected for 
crosstalk in order to stay on the white cut shown in Fig. 2a, 
but the actually measured Vsgi is shown here. Arrows denote 
the transition region between and vr-behavior, associated to 
the crossing point of the switching current. 

The traces exhibit two high switching current peaks 
corresponding to the Coulomb degeneracy points on the 
sides of the Kondo ridge in an odd valley of QDl. Record- 
ing such traces at different magnetic fields provides ac- 
cess to the flux modulation of the switching current in the 
nano- SQUID. Increasing the magnetic flux (p from to 
^o/2 leads to a steady decrease of I^w outside the odd oc- 
cupancy region of JJl, which corresponds to a standard 
0-type behavior^'^ in the Coulomb blockaded even valleys 
of QDl, see Eq. 1. The flux dependence of the switching 
current within the odd-charge Kondo domain turns out 
to be more interesting, as we will analyze now. Clearly, 
the magnetic field behavior of I^w is reversed deep inside 
the odd occupancy region of QDl, as the switching cur- 
rent is greater for cj) = 0o/2 than for = 0, indicating 
a TT-type Josephson behavior. One can therefore iden- 
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FIG. 6: O-TT phase transition diagram and nano-SQUID 
modulations, a) The experimentally determined phase 
boundary between and vr-j unction behavior of JJl (orange 
squares) is given as a function of the dot energy cq/U and the 
level width F/U. The Coulomb repulsion U ^ O.SmeV was 
estimated from the finite bias spectroscopy of the Coulomb 
blockade diamond. Error bars indicate the uncertainty in the 
estimate of F from the analysis of the Kondo resonances. 
The purple line represents the theoretical phase diagram 
for U/A = 10 corresponding to the experimental value of 
A ^ 80/ieV, see Appendix B 2. b) Magnetic field modula- 
tions of the nano-SQUID switching current taken in the mid- 
dle of the odd-charge Kondo ridge (red star in Fig. 2a) for 
different backgate voltages associated to the fine mesh of dots 
at eo = in panel a). For better comparison, the switch- 
ing current modulations are all normalised to the maximum 
current amplitude, which is strongly suppressed in the tran- 
sition region. A clear non-harmonic regime occurs near the 
O-TT phase boundary, where bistable behavior (in phase) of the 
supercurrent may be attributed to the first order transition 
between the and tt states. 



tify precisely from the crossing of the switching current 
traces at which sidegate voltage Vsgi (related to the dot 
energy) the behavior changes from to tt type. This al- 
lows to define a O-tt phase boundary for a given backgate 
voltage. Now, by similarly examining the /gw charac- 
teristics at different backgate voltages (allowing to tune 
the linewidth F), we note that decreasing the backgate 
voltage {i.e. enhancing F) reduces the range for tt be- 
havior, until the tt phase completely collapses below the 
critical Vbg = —0.65V and a 0-junction is maintained 
all along the Kondo ridge. This physical behavior can 
be expected from the stronger Kondo screening at larger 
F that tends to favor the 0-state. From these measure- 
ments, we can unambiguously assign a or tt behav- 
ior to the JJl, as a function of both the level position 
eo and the width F of QDl, as determined previously 
from the analysis of the normal state transport data. 
For all recorded transitions (corresponding to the black 
arrows in Fig. 5), we have extracted the corresponding 
microscopic parameters F and eo, and plotted them on 
an experimental phase diagram shown in Fig. 6a. As a 
quantitative test of our analysis, we have displayed on 
Fig. 6a the theoretical phase diagram obtained from a 
self-consistent description of Andreev bound states^^ for 
/7/A ^ 10 which corresponds to the experimentally mea- 
sured U O.SmeV and A 80/ieV. Error bars are tak- 
ing into account the uncertainty in the determination of 
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r from the finite bias Kondo resonances, see Sec. II C. 
Tiie bell-siiape of tiie piiase boundary togetiier witii tiie 
nearly quantitative comparison to tiieory gives strength 
to the interpretation of the O-tt transition as a first order 
phase transition associated to the crossing of the Andreev 
bound states at the Fermi level^^'^^. A key point here is 
that Kondo screening is decisive to allow the existence 
of the 0-phase in the center of the odd charge Coulomb 
diamond in our experimental conditions {U GF), see 
Fig. 7 in Appendix A. While the O-tt transition is always 
related to a simple Andreev level crossing, this analysis of 
our data clearly demonstrates that it is the competition 
between the normal state Kondo temperature Tk and 
the superconducting gap A which determines the precise 
location of the O-tt phase boundary^^'^^. 

From theoretical expectations^^'^^, a second possible 
smoking gun for the first order O-tt phase transition lies 
in the anharmonic behavior (in phase) of the Joseph- 
son junction in close vicinity to the O-tt phase boundary. 
This prediction motivates us to consider the field mod- 
ulation of the switching current /gw with fine changes of 
the backgate voltage. This is plotted on Fig. 6b, where 
the quantum dot level was taken in the center of the 
odd Kondo valley (particle- hole symmetric point eo = 0, 
corresponding to the red star in Fig. 2a and to the fine 
mesh of dots in Fig. 6a). For Vbg > —0.6V modula- 
tions show that JJl is a 7r-junction because of the (j)o/2 
shift from those of a normal SQUID. On the contrary 
for Vbg < —0.7V, modulations turn back to the stan- 
dard behavior, indicating that JJl is a 0-j unction. How- 
ever, for —0.6 > Vbg > —0.7 the nano-SQUID switching 
current modulations show strong anharmonicities. For 
this range of backgate voltage, the critical current of JJl 
is very small, leading to a strongly asymmetric SQUID, 
thus JJ2 always switches at the same phase difference, 
implying that /gw refiects directly the current-phase re- 
lation of JJl^^. The observed non-harmonic signal could 
be interpreted as a further indication for the bistable be- 
havior of the junction associated with the first order O-tt 
transitioni8-20,23,27_ 



IV. CONCLUSION 

In conclusion, we have realized a nano-SQUID based 
on a superconducting carbon nanotube quantum dots 



that is fully tunable thanks a set of electrostatic gates, 
allowing a precise control of the microscopic parameters 
of the device. This allowed us to determine an experi- 
mental phase diagram for the O-tt transition in the Kondo 
regime, which turned in good agreement with theoreti- 
cal calculations based on the competition between the 
Kondo temperature and the superconducting gap. The 
observation of anharmonic behavior in the supercurrent- 
phase relation near the phase boundary is consistent with 
the first order nature of the O-tt transition associated to 
the crossing of Andreev levels. Fascinating prospects of- 
fered by the present work are the control and monitoring 
of the O-TT transition from supercurrent measurements, as 
performed here, with simultaneous local spectroscopy of 
the Andreev spectrum on the quantum dots in the spirit 
of the recent measurements of Fillet et al.^^. Such future 
developments of our experiment, which seems achievable 
by probing the nano-SQUID with a scanning tunneling 
microscope, would bring strongly correlated supercon- 
ducting nanostructures to a new level of control and un- 
derstanding. 
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Appendix A: Theoretical analysis of the — tt phase 
diagram 



1. Model for a superconducting quantum dot 

The standard Hamiltonian to describe a single super- 
conducting quantum dot is given by the superconducting 
Anderson model 



(Al) 
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where 

Hi = yercl C-, .-V fA,4^y -,, . + h.c.) 

k,a k 

Ha = ^(eo + t^/2)4c«a + t^W 

k^cr 

In the above equations, is the annihilation operator of 
an electron with spin a on the dot, c-* that of an elec- 

tron with spin a and wave vector k in the lead i = L^R^ 
and Ua- = dl^d^. The leads are described by standard 
s-wave BCS Hamiltonians Hi with superconducting gaps 
= A e*'^^ . The phase difference of the latter is noted 
(f = (fL — ^R- Furthermore, the leads are assumed to 
have flat and symmetric conduction bands, i.e. the ki- 
netic energy e- . measured from the Fermi level ranges in 

[—D^D] and the density of states is po = 1/(2D). We 
assume A:- independent and symmetric tunneling ampli- 
tudes t between the dot and both superconducting leads. 
The dot is described by a single energy level eo submit- 
ted to the Coulomb interaction U (in our convention eo 
vanishes at the center of the Coulomb diamond). 

2. Renormalized Andreev Bound states and phase 
diagram of the O-tt transition 

In the superconducting state, the four atomic states 
of the quantum dot evolve onto renormalized Andreev 



bound states (ABS) that possibly live within the gap. A 
quantitative description of this process was proposed in 
Ref. 28, starting with bare values of the ABS splitting in 
the limit of infinite gap: 



SEl = El-El = ^ + ^e,-+T^' . (A3) 



with 



= r- arctan ( ^ ) cos (^^ j . (A4) 



In this simplified (and unrealistic) limit, the O/tt transi- 
tion corresponds to the crossing of the |— ) and |cr) states, 
which occurs for 5E^ = 0, leading to a dome-like shape 
in the {eo/U ,T /U) plane. However, the phase boundary 
quantitatively depends on the precise value of the su- 
perconducting gap A, which must be more realistically 
included in the calculation. This is done by calculating 
the corrections at order 1/A to the ABS positions^^, fol- 
lowed by a self- consistency loop that takes into account 
the leading logarithmic singularities: 



2A 



and 



^ Jo 

(I) 



cos 



5E[ 

2A 

~E 



(I) 



uv 



COS 



E - SE_ (A) E -h SEl E -h SE^ 



E-SE_{A) E^SEl E^SE^ 



_E-SE^{A) E^SEl E^SE^ 
-2 1 1 



E-SE^{A) E^SEl E^SE^ 

I 



■2\T^\uv 



(A5) 



2\r^\uv , 



(A6) 



with E = Ve^TA^, and 5E^_,5E% have been defined 
in Eqs. (A2)-(A3). The numerical resolution of the 
self-consistent equation (A5) provides an accurate de- 
termination of the phase boundary under the condition 
SE-{A) = 0, that we successfully compared to the ex- 
perimental data. Note that we corrected here a misprint 



in Ref. 28, namely a factor 2 in front of the second term 
within the integral in (A5)-(A6). 

In order to stress the key role of the Kondo effect 
for the O-TT transition in our experimental conditions, 
we compare the phase diagram at particle-hole symme- 
try (eo = 0) obtained from the renormalized-ABS the- 
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Renormalized ABS 
Hartree-Fock 
Experimental range 




FIG. 7: Theoretical phase diagram for the O-tt transi- 
tion at particle-hole symmetry. A comparison between 
the renormalized- ABS theory (which includes Kondo correla- 
tions) and static Hartree-Fock theory is made. Computations 
were done for the large bandwidth limit of the superconduct- 
ing Anderson model Eq. (Al), and the experimental range of 
the operating region of Figs. 5 and 6 was added for clarity. 
This comparison shows the key role of Kondo screening to 
allow the existence of a O-tt transition at intermediate corre- 
lations (U > ttF). 



ory^^ and from static Hartree-Fock mean- field theory^^, 
see Fig. 7. Because the renormalized-ABS approach in- 
cludes the Kondo scale (at one- loop order), it allows the 
extension of the O-tt boundary for arbitrary large val- 
ues of U/T. In contrast, the static mean- field approach 
is unable to restore a 0-state for Coulomb interaction 
such that [/ ^ ttF, and fails to reproduce our experimen- 
tal observation of a supercurrent reversal in the regime 
f/ :^ 6r. This comparison shows that the phase bound- 
ary in our experiment is indeed associated to a compe- 
tition between the normal state Kondo temperature and 
the superconducting gap, in agreement with theoretical 
expect at ions^^'^^ . 



Appendix B: Determination of the microscopic 
parameters of the nanosquid 



1. Charging energy 



Because the charging energy in a carbon nanotube 
quantum dot results from the confinement between fixed 
contacts, one does not expect large variations of U for 
small detuning of the backgate. In order to determine 
the experimental phase diagram for the — tt transition, 
an estimate of U is required. This is obtained by con- 
sidering the Coulomb stability diagram of JJl for two 
different values of the backgate, see Fig. 8, and extrap- 
olating the diamond edges to large bias. Because of the 
large linewidth of our strongly coupled nanostructure, 
the determination of U also contains sizeable error bars, 
which we estimate as = 0.80 ± O.OSmeV. 
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FIG. 8: Conductance map of JJl for two values of 
the backgate voltage. Extending the diamond edges to 
finite bias allows to extract the Coulomb repulsion of the dot 
U = 0.80 =b 0.05meV. 



2. Proximity gap 

Current-bias measurements were performed in order to 
directly access both the superconducting switching cur- 
rent and the differential conductance of the nano-SQUID 
at T = 35mK. In the presence of superconductivity, 
the two cotunneling peaks associated with quasiparti- 
cle current in the differential conductance^^ appear at 
V = ±2A/e ~ ±160/iV, where 2 A is the superconduct- 
ing gap provided by the proximity effect on the nanotube, 
see Fig. 9. This allows to extract the superconducting 
gap in our device, A ^ 80/ieV, which is reduced from 
the bulk value of Abuik = 175/ieV for aluminum, due 
to the thin palladium contact layer between the carbon 
nanotube and the aluminum electrodes. 
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FIG. 9: Signature of the superconducting 50 nm thick 
aluminum leads. The differential conductance vs bias volt- 
age shows the presence of a superconducting gap around zero- 
bias, as well as the cotunneling peaks arising from the quasi- 
particule tunneling. The gap value A ^ 80 /xeV is thus ob- 
tained. The measurement was done in a blocked region of the 
nano-SQUID, which explain that no supercurrent is visible at 
zero bias. 



Appendix C: Switching current detection 

In order to have an accurate method to extract the 
switching current from volt age /current characteristics, 
even for small transition voltage jumps to the dissipative 
state, we have implemented a digital filter based on the 
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FIG. 10: Digital filter for switching current detection. 

a) Schematic view of the implemented filter based on the 
work of Liu et al.^^ . The first order moment /i(t) of the sig- 
nal x{t) in a sliding window is obtained through the filter 
hi(t) = Rect(t/Li)/Li, with Rect(t) the normalised rectan- 
gular function and Li the filter length. The variance cr^(t) 

of /i(t) is simply calculated as (/i(t)^)-(/i(t)) in a sliding 
window with /i2(t) = Rect{t/ L2) / L2. b) Voltage/Current 
characteristics and the estimated variance of the first order 
moment obtained by the implemented digital filter. 



work of Liu et al.^^. The main purpose of this filter is the 
detection of transitions from noisy signal, which we apply 
to the superconducting/ normal transition. The opera- 
tion consists in estimating the variance of the first order 
moment of the signal in a sliding window. A schematic 
view of the filter is presented in Fig. 10a. The first or- 
der moment is estimated to begin with via a classic 
averaging filter in a sliding window characterized by the 
impulse response hi{t) = Rect(t/Li)/Li with Rect(t) the 
normalised rectangular function and Li the filter length. 
Finally the estimated variance is obtained by - 

{lJi(t))^ with another averaging filter h2(t) of length L2. 
For the switching current detection with a sample rate 
of one thousand, we have taken Li = L2 = 4. Such filter 
provides a sharp signal from the steplike features of our 
volt age /current characteristics as presented n Fig. 10b. 



